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Abstract 

We show that the density of quadratic forms in n variables over Zp that are isotropic is a rational 
function of p, where the rational function is independent of p, and we determine this rational func¬ 
tion explicitly. When real quadratic forms in n variables are distributed according to the Gaussian 
Orthogonal Ensemble (GOE) of random matrix theory, we determine explicitly the probability that 
a random such real quadratic form is isotropic (i.e., indefinite). 

As a consequence, for each n, we determine an exact expression for the probability that a 
random integral quadratic form in n variables is isotropic (i.e., has a nontrivial zero over Z), when 
these integral quadratic forms are chosen according to the GOE distribution. In particular, we find 
an exact expression for the probability that a random integral quaternary quadratic form is isotropic; 
numerically, this probability of isotropy is approximately 98.3%. 


1 Introduction 

An integral quadratic form Q inn variables is a homogeneous quadratic polynomial 

Q{Xi,X2, ...,Xn)= ^ CijXiXj, (1) 

where all coefficients Cij lie in Z. The quadratic form Q is said to be isotropic if it represents 0, i.e., if 
there exists a nonzero n-tuple (fci,..., e Z” such that Q{ki^... ^ k„) = 0. We wish to consider the 
question: what is the probability that a random integral quadratic form in n variables is isotropic? 

In this paper, we give a complete answer to this question for all n, when integral quadratic 
forms in n variables are chosen according to the Gaussian Orthogonal Ensemble (GOE) of random 
matrix theory [1]. In particular, in the most interesting case n = 4, we show that the probability that a 
random integral quaternary quadratic form is isotropic is given by 

(2 + n “ 4(p + l)2(p4 + p3+p2 + p+l)) ~ ( 2 ) 

More precisely, let D be a piecewise smooth rapidly decaying function on the vector space 
]^n(n+i )/2 qP quadratic forms in n variables (i.e., D{x) and all its partial derivatives are o(|a:|“^) 
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for all N > 0), and assume that D{Q)dQ = 1; we call such a function D a nice distribution on the 
space of real n-ary quadratic forms. Then we define the probability, with respect to the distribution D, 
that a random integral n-ary quadratic form Q has a property P by 


lim 

JV —^OO 


-/Q integral, with property P D{Q/X) 


(3) 


if the limit exists. Let denote the probability with respect to the distribution D that a random 
integral quadratic form in n variables is isotropic. If D = GOE is the distribution on the space of 
n X n symmetric matrices given by -f A^), where each entry of the matrix A is an identical 

and independently distributed Gaussian—i.e., the Gaussian Orthogonal Ensemble—then we use pn ■ = 
Pn^^ to denote the probability, with respect to the GOE distribution, that a random n-ary quadratic 
form over Z is isotropic. 

We wish to explicitly determine the probability that a random n-ary quadratic form over Z, 
with respect to the GOE distribution, is isotropic, i.e., has a nontrivial zero over Z. To accomplish this, 
we first recall the Hasse-Minkowski Theorem, which states that a quadratic form over Z is isotropic if 
and only if it is isotropic over Zp for all p and over M. Eor any distribution D as above, let p^ (p) denote 
the probability that a random integral quadratic form, with respect to the distribution D, is isotropic 
over Zp, and let p^(oo) denote the probability that it is isotropic over M (i.e., is indefinite). Then it 
is not hard to show (for the details, see Section 2) that Pnip) = Pnip) is independent of D, and is 
simply given by the probability that a random n-ary quadratic form over Zp, with respect to the usual 
additive measure on Zp*'"”'"^^^^, is isotropic over Zp. Moreover, we will also show in Section 2 that the 
probability p^{oo) that a random integral quadratic form is isotropic over R is equal to the probability 
that a random real quadratic form (with respect to the same distribution D) is indefinite. 

Eor any distribution D as above, it can be deduced from the work of Poonen and Voloch [6], 
together with the Hasse-Minkowski Theorem, that: 


Theorem 1 The probability p^ that a random (with respect to the distribution D) integral quadratic 
form in n variables is isotropic is given by the product of the local probabilities: 


Pn = Pn{.^)Wpn{p)- 
P 


(4) 


See Section 2 for details. Hence, to determine p^, it suffices to determine Pn(oo) and Pnip) for all p. 

We treat first the probability Pnip) that a random n-ary quadratic form over Zp is isotropic. 
Our main result here is that, for each n, the quantity Pnip) is given by a fixed rational function in p 
that is independent of p (this even includes the case p = 2), and we determine these rational functions 
explicitly. Specifically, we prove the following theorem: 


Theorem 2 Let pnip) denote the probability that a quadratic form in n variables over Zp is isotropic. 
Then 


piip) = 0 , p2ip) = p^ip) = 1 - 2(p^+iy ' 
and Pnip) = f for all n > 5. 


4(p -f l)2(p4 _|_ p3 _|_ p2 _|_ p ’ 
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Our method of proof for Theorem 2 is uniform in n, and relies on establishing certain recursive formu¬ 
lae for densities of local solubility for certain subsets of n-ary quadratic forms defined by their behavior 
modulo powers of p. In particular, we obtain a new recursive proof of the well-known fact that every 
n-ary quadratic form over Qp is isotropic when n > 5. See Section 3 for details. 

We turn next to the probability Pn{oo) = p^°^(oo) that a real n-ary quadratic form is isotropic. 
Closed form expressions for Pn(oo) for n < 3 were first given by Beltran in [4, (7)], while in [2] it was 
argued that 1 — Pn{oo) decays like ('°g3)/4 as n —)■ oo. 

In Section 4, we show how to obtain an exact formula for pn(oo) for any given n. More 
precisely, using the de Bruijn identity [3] for calculating certain determinantal integrals, we express 
Pn(oo) as the Pfaffian of an explicit n' x n' matrix, where n' := 2[n/2], whose entries are given in 
terms of values of the gamma and incomplete beta functions at integers and half-integers. Indeed, let 
r denote the usual gamma function r(s) = x^~^e~^dx, and {3t the usual incomplete beta function 

f^tih j) = Jo — xy~^dx. Then we have the following theorem giving expressions for p„(oo): 


Theorem 3 Let n > 1 be any integer, and define n' := 2[n/2]. When real n-ary quadratic forms 
are chosen according to the n-dimensional Gaussian Orthogonal Ensemble, the probability of isotropy 
over M is given by 


Pn(oo) 


Pf(2l) 


2(n-l)(n+4)/4 r(f) 


(5) 


where A is the n' x n' skew-symmetric matrix whose {i,j)-entry atj is given for i < j by 


2.+,-2r(!M) (i, 1) _ ^.(1, d) n, 

2*-^r(i) if i<j = n + l. 


( 6 ) 


{Note that the second case in (6) arises only when n is odd.) 


Theorem 3 allows one to calculate Pn(oo) exactly in closed form for any given n. In particular, 
it follows from the Pfaffian representation in Theorem 3 that pn(oc) is a polynomial in 7r~^ of degree at 
most with coefficients in Q(\/2) (see Remark 9). In Table 1, we give the resulting formulae for 
Pn(oo) for all n < 8, and also provide numerical approximations. (For any n > 8, we have Pn(oo) ~ 1 
to more than 10 decimal places!) 

Combining Theorems 1, 2, and 3, we finally obtain the following theorem giving the probability 
Pn that a random integral quadratic form in n variables has an integral zero: 


Theorem 4 Let D be any nice (i.e., piecewise smooth and rapidly decaying) distribution. Then the 
probability p^ that a random integral quadratic form in n variables with respect to the distribution D 
is isotropic is given by 


0 


Pn = 


pf(oo)n (i 


po 


ifn < 3; 
ifn = 4; 


4(p -f l)2(p^ + p^ + p^ + p + 1) 

p^(oo) ifn > 5. 

IfD = GOE is the GOE distribution, then the quantities Pn{oo) = p^(cxd) are as given in Theorem 3. 
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n 

Pn(oo) = 

Pn(oo) 

1 

0 

0 

2 

\/2/2 

0.7071067811 

3 

1/2 + x/2 7r-i 

0.9501581580 

4 

1/2 + y/2/8 + 71-^ 

0.9950865814 

5 

3/4+ (2/3 + \/2/12)7r-^ 

0.9997197706 

6 

3/4 + 7\/2/64 + (37/48 - ^2/3)7^^ 

0.9999907596 

7 

7/8 + (47/120 + 109V2/480)7r-^ - (32^2/45)77-2 

0.9999998239 

8 

7/8 + 9\/2/256 + (2377/3840 - 53\/2/480)7r-i - (32/45)7r-2 

0.9999999980 


Table 1: Probability Pn(oo) that a random n-ary quadratic form over M from the GOE distribution is 

isotropic, for n < 8. 


In particular, when B = GOE, we have = 0 for n = 1, 2, and 3, while for n = 4 we obtain 
the expression (2) for p 4 . For n > 5, we have pn = Pn{oo), and so the values of are as given by 
Theorem 3. Theorem 4 shows that n = 4 is in a sense the most interesting case, as all places play a 
nontrivial role in the final answer. 

It is also interesting to compare how the probabilities change if instead of the GOE we use 
the uniform distribution U on quadratic forms, where each coefficient of the quadratic form is chosen 
uniformly in the interval [—1/2,1/2]. While the quantities p]^(oo) can easily be expressed as explicit 
definite integrals, it seems unlikely that they can be evaluated in compact and closed form for general n 
in this case. Using numerical integration, or a Monte Carlo approximation, we can compute /o]^(oo) ~ 
0, 0.627, 0.901, 0.982, 0.998, and > 0.999 for n = 1, 2, 3, 4, 5, and 6, respectively. It is known (see, 
e.g., [I, Theorem. 2.3.5]) that 1 — Pn(oo) decays faster than for some constant c > 0; the actual 
rate of decay is likely even faster. 

In particular, we have p^ = pf (oo) Hp P4(p) ~ 97.0%, which is slightly smaller than the GOE 
probability pf^^ ~ 98.3%. We summarize the values of p^, and provide numerical values in the cases 
of the uniform and GOE distributions, in Table 2. 

This paper is organized as follows. In Section 2, we prove the product formula in Theorem 1. 
The theorem is known in the case of the uniform distribution U (or indeed any uniform distribution 
supported on a box) for any n > 4 by the work of Poonen and Voloch [6], which in turn depends on the 
Ekedahl sieve [5]. To complete the proof of Theorem 1, we first prove directly that both sides of (4) 
are equal to 0 for n < 3. For n > 4, we prove that (4) is true for a general nice distribution B by 
approximating U by a finite weighted average of uniform box distributions, where the result is already 
known. The condition that B is rapidly decreasing (as in the case ofB = GOE) plays a key role in the 
proof; indeed, we show how counterexamples to (4) can be constructed when this condition does not 
hold. 

In Section 3, we then prove Theorem 2, i.e., we determine for each n the exact p-adic density of 
n-ary quadratic forms over Zp that are isotropic. The outline of the proof is as follows. First, we note 
that a quadratic form in n variables defined over Zp can be anisotropic only if its reduction modulo p 
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n 

rn 

pU 

rn 

Pn 

1 

0 

0 

0 

2 

0 

0 

0 

3 

0 

0 

0 

4 

/54(oo)np^l 4(p|l)2(p4|p3^p2|p^4) ^ 

97.0% 

98.3% 

5 


99.8% 

> 99.9% 

> 6 

Pn(oo) 

> 99.9% 

> 99.9% 


Table 2: Probability that a random integral quadratie form in n variables is isotropie, for a general 
distribution D, for the uniform distribution, and for the GOE distribution. 


has either two eonjugate linear faetors over Fp 2 or a repeated linear faetor over Fp. We first eompute the 
probability of eaeh of these eases oeeurring, whieh is elementary. We then determine the probabilities 
of isotropy in eaeh of these two oases by developing oertain reeursive formulae for these probabilities, 
in terms of other suitable quantities, whieh allow us to solve and obtain exaot algebraie expressions 
for these probabilities for eaeh value of n. We note that our general argument shows in partioular that 
quadratic forms in n > 5 variables over Qp are always isotropic, thus yielding a new recursive proof 
of this well-known fact. 

Finally, we prove Theorem 3 in Section 4, i.e., we determine for each n the probability that a 
random real n-ary quadratic form from the GOE distribution is indefinite. We accomplish this by first 
expressing, as a certain determinantal integral, the probability that an n x n symmetric matrix from 
the GOE distribution has all positive eigenvalues. We then show how this determinantal integral can 
be evaluated using the de Bruijn identity [3], allowing us to obtain an expression for the probability 
of positive definiteness in terms of the Pfaffian of an explicit skew-symmetric matrix A, as given in 
Theorem 3. We note that the values of these probabilities were known previously for n < 3 (cf. [4, 
(7)]). 

We end this introduction by remarking that the analogues of Theorems 2 and 4 also hold over 
a general local or global field, respectively. Here, we define global densities of quadrics as in [6, 
§4]; more general densities with respect to “nice distributions” could also be defined in an analogous 
manner. Indeed, the analogue of Theorem 1 holds (with the identical proof), where the product on the 
right hand side of (4) should be taken over all finite and infinite places of the number field (the densities 
at the complex places are all equal to 1, since all quadratic forms over C are isotropic). Theorem 2 also 
holds over any finite extension of Qp, with the same proof, provided that when making substitutions in 
the proofs we replace p by a uniformiser, and when computing probabilities we replace p by the order 
of the residue field. 
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2 The local product formula: Proof of Theorem 1 


Let D be any niee (pieeewise smooth and rapidly deeaying) distribution. Our aim in this section is to 
prove the following three assertions from the introduction: 

(a) Pnip) is equal to the probability Pn{p) that a random n-ary quadratic form over Zp, with respect 

to the usual additive measure on is isotropic over Z^; 

(b) p^(cx)) is equal to the probability that a random n-ary quadratic form over M, with respect to the 
distribution D, is indefinite; and 

(c) p^ = pn{oo) ]\ppn{p) (i-e., Theorem 1 holds). 

Items (a) and (b) are trivial in the case that D = \J is the uniform distribution, or more generally 
when D is any distribution U(a, b) that is constant on a box [a, 6] := [oi, 6i] x • • • x [an{n+i)/ 2 , &n(n+i)/ 2 ] 
and 0 outside this box; here a = (oi,..., an{n+i)/ 2 ) and b = (6i,..., &„(n+i)/ 2 ) are vectors in 
such that tti < bi for all i. 

Meanwhile, Theorem 1 for n > 4, in the case that D is the uniform distribution U, follows 
from the work of Poonen and Voloch [6, Theorem 3.6] (which establishes the product formula for the 
probability that an integral quadratic form with respect to the distribution D is locally soluble), together 
with the Hasse-Minkowski Theorem (which states that a quadratic form is isotropic if and only if it is 
locally soluble). In fact, the proof of [6, Theorem 3.6] (which in turn relies on the Ekedahl’s sieve [5]) 
immediately adapts to the case where D = \J{a, b) without essential change. 

To show that Theorem 1 holds also when D = U{a, b) and n < 3, it suffices to prove that in 
this case both sides of (4) are equal to 0. To see this, we may use Theorem 2, which does not rely on 
the results of this section, and which states that the probability that a random n-ary quadratic form over 
Zp is isotropic is equal to Pnip) = 0, 1/2, or 1 — p/(2(p -f 1)^) for n = 1, 2, or 3, respectively. This 
immediately implies that the right hand side of (4) is zero. To see that the left hand side of (4) is zero, 
we note that if a quadratic form over Z is isotropic, then it must be isotropic over Zp for all p (the easy 
direction of the Hasse-Minkowski Theorem). By the Chinese Remainder Theorem, the (limsup of the) 
probability p^ that a random integral n-ary quadratic form is isotropic with respect to the distribution 
H = U(a, b) is at most 

n 

p<Y 

for any F > 0. Letting Y now tend to infinity shows that = 0 for n = 1, 2, or 3, i.e., the left hand 
side of (4) is also zero. 

Thus we have established items (a)-(c), for all n, in the case that D = \J{a, b) is a constant 
distribution supported on a box [a, b]. Clearly (a)-(c) then must hold also for any finite weighted 
average of such box distributions U(-, •). 

To show that (a)-(c) hold for general nice distributions D, we make use of the following ele¬ 
mentary lemma regarding integration of rapidly decaying functions. 
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Lemma 5 Let f be any piecewise smooth rapidly decaying function on Then 


[ fiy)dy= lim^^ 


( 7 ) 


Proof: For any N > 0, let fN^y) be equal to f{y) if \y\ < N, and 0 otherwise. Then /jv is pieeewise 
smooth with bounded support, and so is Riemann integrable. Thus we have 

/ My/^)- (S) 

Sinee / is rapidly deereasing, for any e > 0 we may ehoose N large enough so that \f{y)\dy<e 
and (1/X™) \y/x\>N \fiy/^)\ < ^ > 1- For this value of N, the left hand side of 

(8) is within e of the left hand side of (7), while for eaeh X > 1, the expression in the limit on the right 
hand side of (8) is within e of the expression in the limit on the right hand side of (7). Sinee we have 
equality in (8), we eonelude that the left hand side of (7) is within 2e of both the liminfx^oo and the 
limsupjf^oo of the expression in the limit of the right hand side of (7). Sinee e is arbitrarily small, we 
have proven (7). □ 

Note that Lemma 5 does not neeessarily hold if we drop the eondition that / is rapidly deeay- 
ing. For example, if / the eharaeteristie of a finite-volume region having a eusp going off to infinity 
eontaining a rational line through the origin (and thus infinitely many lattiee points on that line), then 
the left hand side of (7) is finite while the expression in the limit on the right hand side of (7) is infinite 
for any rational value of X. 

Lemma 5 implies in partieular that 

^ = l (9) 

Q integral 


for any niee distribution D. 

Now any pieeewise smooth rapidly deeaying funetion ean be approximated arbitrarily well by 
a finite linear eombination of eharaeteristie funetions of boxes. Let D be a niee distribution. For any 
£ > 0, we may find a niee distribution that is a finite weighted average of box distributions U(-, •), 
sueh that 

J\D{y) - D,{y)\dy < e. (10) 

By Lemma 5, we then have 

Xn{n+l)/2 E \D{Q/X)-D,{Q/X)\<e. (11) 

Q integral 
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To show that p^(p) = pn{p), we note that 


Pn{p) = Pn"iP) = lim 

JV^CXD 

XIq integral,isotropic/Zp 

(12) 

Eein..g,.,£>.((?/ a:) 

integral,isotropic/Zp ^£{.Q / ^ ) 

= lim 

X — ^ oo 

(13) 

/^n{n+l)/2 

= lim 

EQintegral,isotropic/Zp X) + ^Xi 

(14) 

J^n(n+l)/2 

X —^oo 

= lim 

EQintegral,isotropic/Zp ^{QX) + ^X^ 

(15) 


X — ^OO 


where for suffieiently large X we have \E{X, e) \ < by (11); here the first equality follows 

beeause is a finite weighted average of box distributions U(-, •), the second equality follows from 
the definition (3), and the third and fifth equalities follow from (9). Letting e tend to 0 in (15) now 
yields pnip) = Pnip)^ proving item (a) for general nice distributions D. 

Analogously, we have 


' Q isotropic/R 


De{Q)dQ = p, 


D, 


[oo) = lim 

X^oo 


= lim 

X^OO 


integral,isotropic/R DAQ/X) 

integral,isotropic/M-^(^/^) + ^ Xx) 
Eg integral ^ 


(16) 

(17) 


where again for sufficiently large X we have \E'{X, e) \ < By (10), the leftmost expression 

in (16) approaches Jq isotropic/R E(Q)dQ as e —)■ 0, while expression (17) approaches p^(cxo) by defi¬ 
nition (3). This thus proves item (b) for general nice distributions. In particular, we have also proven 
that 

limp/"(«)) = p/(oo). (18) 

£^•0 

Finally, we have in a similar manner: 


Pn 


/ \ TT / \ 1 - 5^0 integral,isotropic/Z 

(oo) n p-(p)=ft. = i- 


= lim 

X —^OO 


EQintegral,isotropic/Z ^{QX) + ^ Xi 


E 


Q integral 


D{QIX) 


(19) 

( 20 ) 


where again for sufficiently large X we have \E''{X, e) | < _ gy (18), the leftmost expression 

in (19) approaches PnX) rip Pn{p) as e —)■ 0, while expression (20) approaches p^ by definition. We 
have proven also item (c) for general nice distributions, as desired. 
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3 The density of n-ary quadratic forms over that are isotropic 
Proof of Theorem 2 


3.1 Preliminaries on n-ary quadratic forms over 

Fix a prime p. For any free Zp-module V of finite rank, there is a unique additive p-adic Haar measure 
Pv on ^ whieh we always normali z e so that pv{y) = 1- All densities/probabilities are eomputed with 
respeet to this measure. In this seetion, we take V = Vn to be the n{n + 1)/2-dimensional Zp-module 
of n-ary quadratie forms over Zp. We then work out the density Pn{p) (i-C- measure with respeet to pv) 
of the set of n-ary quadratic forms over Zp that are isotropic. 

We start by observing that a primitive n-ary quadratic form over Zp can be anisotropic only 
if, either: (I) the reduction modulo p factors into two conjugate linear factors defined over a quadratic 
extension of Fp, or (II) the reduction modulo p is a constant times the square of a linear form over Fp. 
Let and Q be the probabilities of Cases I and II, i.e. the densities of these two types of quadratic 
forms in 14- Then 

An) _ T _ An) _ An) _ ^ 

'^0 SI S2 pn{n+l)/2 

is the probability that a form is primitive, but not in Cases I or II. Let (resp. be the probability 
of isotropy for quadratic forms in Case I (resp. Case II). Then 


in) I AA An) 


Pn{p)=ir+i\ 


a 


1 


+ A A + 


pn{n+l)/2 


Pn(p), 


implying that 


Pn{p) = 


P' 


n(n+i)/2 


P 


,n{n+l)/‘2 


(A + AA+AA)- 


in)An) I A'^AAP 


( 21 ) 


3.2 Some counting over finite fields 


Let (resp. be the probability that a quadratic form is in Case I (resp. Case II) given the “point 
condition” that the coefficient of xl is a unit. Similarly, let be the probability that a quadratic form 
is in Case I given the “line condition” that the binary quadratic form Q{xi,X 2 , 0,..., 0) is irreducible 
modulo p. Note that it is impossible to be in Case II given the line condition, but we may also define 

The values of 
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A) 


= 0. Set = 1 — pj"'' — pg"''' and z/q"'' = 1 — z/j"'' — = 1 — u 

, are given by the following easy lemma. 




An) 


(n) 


An) 


An) 


An) 


Lemma 6 The probabilities that a random quadratic form over Zp is in Case I or Case II are as 
follows. 

• Case I (all; relative to point condition; relative to line condition) 

An) ^ (p’"-l)(p^-p) _ (n) ^ P'^~^ - 1 . (n) ^ 1 

2(p-f 2pnA-PA' p{n-l){n-2)/2' 
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• Case II (all; relative to point condition; relative to line condition) 


An) 

S2 


_ 

pn{n+l)/2 ’ 



1 

pn(n—l )/2 ’ 



= 0 . 


Proof: Case I: There are — l)/(p^ — 1) linear forms over Fp 2 up to sealing; subtraeting the 
(pn _ whieh are defined over Fp, dividing by 2 to aceount for eonjugate pairs and then 

multiplying by p — 1 for sealing gives 2 (p+^i)~^^ Case I forms, and henee the value of 

Similarly, the number of Case I quadratie forms satisfying the point condition is (pAn-^) _ 
^n- 1 ) (p _ 1 ) /2 Dividing by the probability 1 — 1 /p of the point condition holding gives p" (p”“^ ~ 1) /2 
and hence the value of p} . 

Lastly, the number of Case I quadratic forms satisfying the line condition is p^”“^(p — 1)^/2; 
dividing by the probability of the line condition holding gives and hence the value of 

Case II is similar and easier: the number of Case II quadratic forms is p"^ — 1, of which p^—pn-i 
satisfy the point condition and none satisfy the line condition; the given formulae follow. □ 


3.3 Recursive formulae 

We now outline our strategy for computing the densities Pn{p) using (21), by evaluating for j = 
1,2. If a quadratic form is in Case I, then we may make a (density-preserving) change of variables, 
transforming it so that its reduction is an irreducible binary form in only two variables. Now isotropy 
forces the values of those variables, in any primitive vector giving a zero, to be multiples of p; so we 
may scale those variables by p and divide the form by p. Similarly, if a form is in Case II, then we 
transform it so that its reduction is the square of a single variable, scale that variable and divide out. 
After carrying out this process once, we again divide into cases and repeat the procedure, which leads 
us back to an earlier situation but with either the line or point conditions, which we need to allow for. 

To make this precise, we introduce some extra notation for the probability of isotropy for 
quadratic forms which are in Case I or Case II after the initial transformation: let (resp. 
be the probability of isotropy given we are in Case I (resp. Case II) after one step when the original 
quadratic form was in Case I, and similarly 7) (resp. 7) ) the probability of isotropy given we are in 
Case I (resp. Case II) after one step when the original quadratic form was in Case II. 


Lemma 7 . 


1 . 


a 


( 2 ) 


= 0, and for n >3, 


(n) _ f(n-2) 


a^ = 


so 






+ 




(n-l)(n-2)/2 


pi 


iAAn)^ 


2. a 2 ^ = 0, and for n > 2, 


4”’ = «L‘’+?C‘’7l”’+eL‘'7f’ 


pn(n—l)/2 
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(2) 

Proof: We have = 0 sinee a binary quadratie form that is irredueible over Fp is anisotropie. Now 
assume that n > 3, and (for Case I) Q{xi ,..., Xn) (mod p) has two eonjugate linear faetors. Without 
loss of generality, the reduetion modulo p is a binary quadratie form in xi and X 2 . Now any primitive 
veetor giving a zero of Q must have its first two eoordinates divisible by p, so replaee Q{xi,, Xn) 
by ^Q{pxi,px 2 , Xs,..., Xn). The reduetion modulo p is now a quadratie form in xs,..., If the 
new Q is identieally zero modulo p, then, after dividing it by p, we obtain a new integral form that 
lands in Cases I and II with probabilities and respeetively, sinee it satisfies the line eondition; 
otherwise, we divide into oases as before, with the probabilities of being in eaoh ease given by 

The result for is proved similarly: without loss of generality the reduetion modulo p is a 
quadratie form in xi only, we replaee Q{xi ,..., x„) by ^Q{pxi,X 2 ,..., x„), whose reduetion mod¬ 
ulo p is a quadratie form in X 2 ,..., Xn. If the new Q is identieally zero modulo p, then, after dividing 
by p, we have an integral form that lands in Cases I and II with probabilities pi and > respeetively, 
sinee it satisfies the point eondition; otherwise, we divide into oases, with probabilities ■ □ 

It remains to oompute (for n > 4), (5^"^ (for n > 3), (for n > 3) and 72 ”^ (for n > 2). 
Sinee = 0, we do not need to oompute or whioh are in any ease undefined. 


Lemma 8 . 

(i) Ifn>4: then -f z/j" also, = 0. 

(ii) Ifn > 3 then (3^'^ = also, (3^'^ = 0. 

(iii) Ifn > 3 then also, y^^ = 0. 

(iv) Ifn > 2 then y^”^ = also, y® = 0. 


Proof: In Case I, the initial transformation leads to a quadratie form for whioh the valuations of the 
ooeffioients satisfy* 


> 

1 

> 

1 

> 

1 . 

. > 

1 

> 

1 

> 

1 

> 

1 . 

. > 

1 

> 

0 

> 

0 

> 

0 . 

. > 

0 



> 

0 

> 

0 . 

. > 

0 





> 

0 . 

. > 

0 


( 22 ) 


> 0 


and (resp. are the probabilities of isotropy given that the reduetion modulo p of the form in 
xs, X 4 ,..., Xn is in Case I (resp. Case II). 

*In this and the similar arrays which follow, we put into position {i,j) the known condition on v(aij), so the top left 
entry refers to the coefficient of xf the top right to xiXn and the bottom right to x^- 


II 






Similarly, in Case II the initial transformation leads to 


> 

1 

> 

1 

> 

1 . 

. > 

1 

> 

0 

> 

0 

> 

0 . 

. > 

0 

> 

0 

> 

0 

> 

0 . 

. > 

0 



> 

0 

> 

0 . 

. > 

0 





> 

0 . 

. > 

0 


> 0 


(23) 


and 7 ^^”^ (resp. are the probabilities of isotropy given that the reduction modulo p of the form in 
0 : 2 , X 3 ,..., is in Case I (resp. Case II). 

(i) To evaluate (3^^ we may assume, after a second linear change of variables, that we have 


> 1 

> 1 

> 1 

>1 ... 

> 1 

> 1 

> 1 

> 1 

>1 ... 

> 1 


> 0 

> 0 

>1 ... 

> 1 



> 0 

>1 ... 

> 1 




>1 ... 

> 1 


> 1 


and that the reductions modulo p of both ^Q(xi, X 2 , 0 ,..., 0 ) and (^( 0 , 0 , 0 : 3 , 0 : 4 , 0 ,..., 0 ) are irre¬ 
ducible binary quadratic forms. Any zero of Q must satisfy X 3 = = 0 (mod p). This gives a 

contradiction when n = 4, so that Q(a:i,..., X 4 ) is anisotropic, and = 0. Otherwise, replacing 
Q{xi,... ,Xn) by ^Q{xs, X 4 , pxi, px 2 , x^,... ,x„) brings us back to the situation in (22). Now, how¬ 
ever, the line condition holds, so that Cases I and II occur with probabilities and 1 ^ 2 "^”^^ = 0 

instead of and . 


(ii) To evaluate we may assume that the valuations of the coefficients satisfy 


> 1 

> 1 

> 1 

... > 1 

> 1 

> 1 

> 1 

... > 1 


= 0 

> 1 

... > 1 



> 1 

... > 1 




> 1 


and that the reduction modulo p of ^Q(xi,X 2 , 0,..., 0) is an irreducible binary quadratic form. Ifn = 3 

then Q is anisotropic, and = 0. Otherwise, replacing Q{xi, ..., Xn) by ^Q(x 2 , X3,pxi, X4, ..., x^) 
brings us back to the situation in (23) but with the line condition, so that Cases I and II occur with 
probabilities and 1 ^ 2 ””^^ instead of ^ 2 ””^^- 
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(iii) For 7!”^ 


, we may assume that the valuations of the eoeffieients satisfy 


= 1 >1 >1 > 1 
>0 >0 > 1 
>0 > 1 
> 1 


> 1 
> 1 
> 1 
> 1 


> 1 


and the reduetion of Q(0, X 2 , 0 : 3 , 0,, 0) modulo p is irredueible. Any zero of Q now satisfies X 2 = 
X 3 = 0 (mod p). When n = 3 this gives a eontradietion, so Q{xi,X 2 , x^) is anisotropie, and yf ^ = 0. 
Otherwise, replaeing Q(xi,... ,Xn) by ^Q{x 3 ,pxi,px 2 , X 4 ,... ,x„) brings us baek to the situation 

in (22) but with the point eondition, so that Cases I and II oeeur with probabilities and 
(iv) Lastly, for 7 ^”^ we may assume that the valuations of the eoeffieients satisfy 

= 1 >1 >1 • • • >1 

= 0 >1 ... >1 

>1 ... >1 


> 1 - 

If n = 2 then Q{xi,X 2 ) is anisotropie , and 72 ^^ = 0. Otherwise, replaeing Q{xi,... ,Xn) by 
^Q{x 2 ,pxi, X 3 ,..., Xn) brings us baek to the situation in (23) but with the point eondition. □ 


3.4 Conclusion 

Using Lemmas 6 and 8 we ean eompute and 7 ^”^ for j = 1,2 and all n: we first determine /3i from 

Lemma 8 (i), then (3^'^ and 7 !”^^ together using Lemma 8 (ii,iii), and finally 72 ”^ using Lemma 8 (iv). 
The following table gives the result: 



a'"’ 

4"> 

(n) 

7i 

(n) 

72 

n = 2 

— 

— 

— 

0 

n = 3 

— 

0 

0 

1/2 

n = A 

0 

( 2 p+l)/( 2 p + 2 ) 

{p + 2 )/{ 2 p + 2 ) 

1 - (p/(4(p2 + P+ 1)) 

n > 5 

1 

1 

1 

1 


Now, using Lemma 7, we eompute and 



(n) 

a) 

(n) 

a\ 

n = 2 

0 

l/{2p + 2) 

n = 3 

1/ {p + 1) 

{p + 2)/{2p + 2) 

n = 4 

1 - {p^ / {2{p + l){f + p + 1))) 

1 - (p^/(4(p + l)(p^ + p‘^ + P+ 1))) 

n > 5 

1 

1 
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Finally, we compute pn{p) using (21), yielding the values stated in Theorem 2. 

Note that our proof of Theorem 2 also yields a (recursive) algorithm to determine whether a 
quadratic form over Qp is isotropic. Tracing through the algorithm, we see that, for a quadratic form of 
nonzero discriminant, only finitely many recursive iterations are possible (since we may organize the 
algorithm so that at each such iteration the discriminant valuation is reduced), i.e., the algorithm always 
terminates. In particular, when n > 5, our algorithm always yields a zero for any n-ary quadratic form 
of nonzero discriminant; hence every nondegenerate quadratic form in n > 5 variables is isotropic. 


4 The density of n-ary quadratic forms over M that are indefinite: 
Proof of Theorem 3 

4.1 Preliminaries on the Gaussian Orthogonal Ensemble (GOE) 

We wish to calculate the probability pn{oo) that a real symmetric matrix M from the n-dimensional 
GOE has an indefinite spectrum. The distribution of matrix entries in the GOE is invariant under or¬ 
thogonal transformations. Since real symmetric matrices can be diagonalised by an orthogonal trans¬ 
formation, the GOE measure can be written directly in terms of the eigenvalues A(M), yielding the 
distribution 

P(A(M) e [A + dA)) = _g|A(A)ine-j"'dA.; (24) 

n 

here 

A(A):= H (A,-A,) = det(yp,(A,)), 

where ((/?j(Aj)) = (A*“^) is a Vandermonde matrix, and the normalizing factor is given by 

z™® = ff Ed. (25) 


See, for example, [1, (2.5.4)]. 

Note that the probability that the matrix M is indefinite is related to the probability that all 
its eigenvalues are positive by 

Pn{oo) = 1 — P(positive definite) — P(negative definite) = 1 — (26) 

where the second equality follows by symmetry. Below we will calculate p^, and hence obtain the 
value of p„(oo). 

4.2 de Bruijn’s identity 

We recall a useful result from [3, §4] for calculating determinantal integrals of the type we will need. 
As a generalisation of an expression for the volume of the space of symmetric unitary matrices, de 
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Bruijn considered integrals of the form: 


n 


det (ifii(xi))dxi... dxr 


a<Xi<-"<Xn<b 


Reeall that the Pfaffian of a skew-symmetrie matrix A = (a^) is given by 

T 


(27) 


(28) 


where r ranges over all partitions 


^ = (*2,72),. • .( 4 , 7 .)} 

of n = 2s where 4 < 4+i and 4 < jk- The sign is of the corresponding permutation 


r = 


12 3 4 

4 7i 4 72 


2 s 

js 


The integral (27) may be rewritten as the Pfaffian of either an n x n skew-symmetrie matrix if 
n is even, or an (n -f 1) x (n -f 1) skew-symmetrie matrix if n is odd. More preeisely, let n' := 2 [?7,/2]; 
then we have = Pf (4), where A is the n' x n' skew-symmetrie matrix whose (i, j)-entry aij is given 
for i < j by 

/a la sign(2/ - x)ipi{x)ipj{y)dxdy if i < j < n; 

r -r • ■ , 1 (29) 

Ja if I < J = n + 1. 


dij 


The second case occurs only when n is odd. Note that this holds for a general measure dx; below, we 
will use dx = e“^^/'^dA, where dA is the Lebesgue measure on M. 

The Pfaffian form of the integral is found by expanding the determinant and using a signature 
funetion to keep traek of the signs and the ordering of the Xi. This signature funetion of n variables ean 
be broken up into a sum of produets of two-variable pieees (and a one-variable piece if n is odd) and 
thus the integral ean be faetorised into a sum of products of two (and one) dimensional integrals whieh 
is reeognised as of the form (28) for a matrix with entries (29). 


4.3 Calculation of p] 


,00 

n 


For a matrix M from the GOE, the joint distribution of the eigenvalues Ai(M) < \ 2 {M) < ■ ■ ■ < 
XniM) is given by 

I 

rlAi<A2<-<A„|A(A)| e 44(1^._ (30) 


^■GOE ■ 


2=1 
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The ordering in the domain of integration below means that we ean replaee |A(A)| by A (A). It then 
follows that is given by the integral 


P 


+ 

n 


^ f...j A(A)n.-WdA. 

"■ 0<Ai<---<An<OO 

J"J det((^i(Aj))]^e-^^»dAi 

"■ 0<Ai<---<An<oo 


(31) 


(*OQ nOQ 


where the last equality follows from the result of §4.2. Here, A = (aij), where for i < j < n we define 

sign(2/ — x)x'^~^y^~^e ^ dxdy 

^ (32) 


Gij — 


JO JO 


i+j 


(H) 


and for n odd we also set ai^n+i = 2*“^r(|). Here the gamma and ineomplete beta funetions are as 
defined in §1. From the resulting skew-symmetrie matrix A, we may evaluate (31) to determine pn(oo), 
yielding Theorem 3. Explieit values of pn(oo) are displayed in Table 1 for n < 8 . 


Remark 9 It is easily shown that the matrix entries in Theorem 3 are of the form x or Xy/ir for 
X G Q(\/ 2 ), in accordanee with whether i + j is even or odd. Let s = [n/2], so that A is a 2s x 2s 
matrix. Then after re-ordering the rows and eolumns we have 


Pf(Al) 


±Pf 




±7r*/2pf 



A 2 

tt-^As 


where Ai, A 2 and A 3 are s x s matriees with entries in Q(\/2). Since 11^=1 r(m/2) = for some 
y G Q, it follows by Theorem 3 and the definition of the Pfaffian that pn(oo) is a polynomial in 
having coefficients in Q(\/ 2 ) and degree at most Ls/ 2 J = |_(n + 1)/4J. 
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